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ABSTRACT 

An electrohydrodynamic traveling wave induction interaction with 
the bulk of a slightly conducting liquid (conductivity of approximately 
10 ^^mhos/meter ) is shovn to produce a fluid flow. A gradient in con- 
ductivity normal to the direction of flow is reqtiired. In the scheme 
described here the conductivity gradient is provided by taking advantage 
of the temperature dependence of conductivity. A traveling potential 
wave is created on an arrangement of electrodes parallel to the flow 
direction and in contact with the liquid. The restilting fields induce 
charges in the bulk which lag the traveling potential wave. Therefore, 
a time -average electric traction in the bulk is created, motivating the 
liquid. Expressions for the fields, the time-average traction, and the 
fliiid velocity are derived and discussed. Experimental results are shown 
which lend credence to the derived equations. 
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CHAPTER I 



INTRODUCTION 

Induced charges in a slightly conducting liquid can interact vith 

traveling wave electric fields to produce a fluid flow. One method of 

producing this kind of electroconvection is to introduce a potential 

wave traveling parallel to an air liquid interface. Induced charges 

in the liquid will relax to the interface and form a traveling wave 

of surface charge which will lag behind the potential wave. An electric 

* 

surface shear will result and the liquid will move. In this kind of 
device the interface is a singular point in a gradient of conductivity. 

An extension of the concept of the surface interaction would be 
to produce a continuous gradient in conductivity in a liquid. Then 
eliminate the interface by bringing the electrodes carrying the travel- 
ing potential wave into contact with the liqiiid. Charges induced in 
the liquid will no longer relax to the surface but will remain in the 
bulk of the liquid. A kind of internal wave of induced charge will 
then travel along behind the traveling potential wave and the possibility 
for an electric traction in the bulk of the liquid is created. 

This phenomenon is very dependent on the charge relaxation time 
which is the ratio of permittivity to conductivity. If the charge 

Melcher, J.R», "Traveling Wave Induced Electroconvection" to be published 
in Physics of Fluids . 
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relaxation time for the liquid is very short, then charge will be 
induced in the liquid but will relax so quickly that it will not 

lag behind the traveling potential wave. On the other hand, if the 
relaxation time is very long, then very little charge will be in- 
duced in the material and the interaction with the . traveling field will 
be negligible. 

EHD induction puit^ing would be feasible for a wide class of 
liqiiids whose relaxation times would be compatible with reasonable 
frequencies for which traveling electric potential waves might be 
generated. 

In succeeding chapters, a theory and experiment will be described 
which deal with the phenomenon of electrohydrodynamic pumping of liqxjids 
in the bulk. The results obtained are somewhat surprising. 
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CHAPTER II 



THEORY 



A. Introduction 

The basis of the following theory for electrohydrodynaraic induction 
pumping in the bulk of a slightly conducting liquid is that the movement 
of charge in the liquid may be explained by a pure conduction model. This 
model applied to induction pumping in the bulk results in currents which 
are small enough that the problem may be considered to be only an electric 
problem and magnetic effects may be ignored. The theory as developed 
here leads to a conclusion which at first glance is somewhat startling 
and violates an intuitive conception of how an induction device ought to 
operate . 

B. Electrical Theory 

1. Electri c Field Equations 

Considering magnetic effects to be small we may immediately write 
Maxwell's equations in the bulk of the liquid. They are: 

V X E = 0 (l) 

^ • eE = q (2) 

* 

Fano, Chu, Adler, Electromagnetic Fields, Energy, and Forces, p.l79« 
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where q is the free charge, s is the permittivity and E is the electric 
field vector. To satisfy conservation of charge 

( 3 ) 

The free current J results from Ohm's Law conduction which written in 
the coordinate system fixed to the laboratory will be 

J = oE + qv (4) 

where a is the electrical conductivity and v is the velocity in the 
fluid flow. Since the Electric Reynolds number given by the formula 

-2 

'Re = cv/cd is small for the problem \inder consideration ('IRe as 2,6 x 10 ) 

the qv term in Equ, (4) may be ignored leading to 

J = oE (5) 

Equation (l) allows that S’ may be written as the negative gradient of 
the potential $. 

E = - ^ $ (6) 

Combining Equs, (2), (3) and (6) leads to: 

9 - + ( V • ( e VO ) - 0 (T) 
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If both cr and e are functions of y in the coordinate system shovn 
in Fig. 1^ this equation is in general difficult to solve^ however^ 
there are sitixations in vhich it may be easily approached. 

One such case occurs if it is assumed that e is a constant 
throughout the bulk of the liquid and that a may be represented as: 

where C is the normalized depth in the channel. C = y/d . re- 
presents some appropriate depth in the channel where a value of a 
is measured which is in some sense an average value of conductivity 
in the channel. It is also necessary that the values of and 
be such that it is reasonable to say that a <j^ in the liquid but 
Vct ~ — cr, /d, where — is the unit vector in the y direction. 

iy iy 

Combining eqviations (8) and (j) with e a constant leads then to 

^o ^ ^ V d^ • ^ ^ ® ^ ^ ^ ° 

5t 

This equation will new, be solved for the two-dimensional case. 

2. Traveling Wave Solutions 

If the upper surface of the channel is excited with a traveling 
potential wave of the form: 

A j(<^ - kx) 

V = Re V^e (lO) 
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Figure 1 Coordinate system in the pumping channel 
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where V is the inagni tude of the applied voltage, vi is the angular 
frequency, k is the wave number, and circumflex , indicates a complex 
qiiantity; then the assumed form of the solution of Equ. (9) is: 

A j(u7fc - kx) 

$ = Re $ (C) e (11) 

Substituting Equ. (ll) into (9) gives: 



a 



A 

i? - 






ac 



(dk) 5 = 0 (12) 



This equation may be solved by assuming for § solution of the form 
. Substituting this solution into Equ. (l2 ) and solving for 
the roots gives: 



a. 



_ 1. 
P - 2 



cr^ + juJe 




o'q + 



f + k(d^f 



(13) 



2 

Since in the experiment performed in Chapter IV, 4(dk) « 

2 

jcue)) this square root may be expanded around the first term. At 
the same time defining new constants T) = ® ~ 

roots may be rewritten; 




_Q_ + ( t + 4(dk)^(i + jS ) ) 

1 + js “ 1 + js n 



(li^) 
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Nov the complete solution for $ may he written and the boundary 
condition on the electric field in the channel may be applied to 
evaluate the constants of integration: 

A P, C PpC 

$ (C) = A^e ^ + A^e ^ (15) 

At the perfectly conducting bottom of the channel where Q = 0, 
the appropriate boundary condition is that the tangential com- 
ponent of electric field must be zero. Therefore, -3$/Sx = 0 
which says that A^ = -A^. At the top of the channel where C = 1 
the condition is that § = V. This condition says that; 




The final form of $ is: 



A P, C PpC 
Re V (e - e ) j(u)t 

§ 2 e 

^1 ^2 
e - e 



- kx) 



(16) 



Since E = 

X 

form; 



S§/5k and = - d$/ddC. , writing E^ and E^ in the 



E = Re E e 

X X 



- kx) 



: E = Re E 6 

y y 



o(u;t - kx) 
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leads to: 



E 

X 



E 

y 



*■ 1*1 ^ Po^ 

ok V (e ■ e ) 

o 

p. Pp 

(e ^ - e ^ ) 

A Pi S vS 
-v„( p,e ^ - p^e M 

a p, p, 

(e ^ - e 



(17) 



(18) 



Equs. ( 17 ) and ( 18 ) complete the solution of the electric field 
problem for the fields in the bulk of the liquid. 

C. Fluid Theory 

1. Force Equation 

The force equation will be written for a two-dimensional 
flow in a re-entrant channel \inder the influence of an internal 
electric shear force which is constant in the x-direction because 
only the time-average value of the electric stress is considered. 
The liquid is incompressible and the flow is steady. 

The general force euqation is : 

p + ( 7 ■ ^ ) 7 - ) + . 7 • i - pg ^ (19) 

5t ^ 



P is the pressure, v the flow velocity, p the density, and g is 
the gravitational acceleration. T. . is the total stress tensor 
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